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ABSTRACT 


A  method  employing  two-dimensional  spline  fits  of  spatially  and  temporally  random  hydro- 
graphic  data  is  developed  in  order  to  be  able  to  determine  seasonally  averaged  geostrophic  cur¬ 
rents.  The  method  is  used  in  an  analysis  of  the  currents  in  the  Southern  California  Bight. 
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THE  DETERMINATION  OF  AVERAGE  GEOSTROPHIC  CURRENT  VELOCITIES  FROM 
TEMPORALLY  AND  SPATIALLY  RANDOM  HYDROGRAPHIC  DATA  WITH  AN  APPLICATION 

TO  THE  SOUTHERN  CALIFORNIA  BIGHT 

INTRODUCTION 

Dynamic  depth  anomaly  data  from  oceanographic  stations  are  available  as  a  standard  product 
from  the  National  Oceanographic  Data  Center  (NODC).  NODC  calculates  these  data  from  tem¬ 
perature  and  salinity  with  depth  data  which  it  obtains  from  various  sources.  It  is  highly  desirous 
to  analyze  these  data  in  a  manner  such  as  to  obtain  seasonally  averaged  dynamic  depth  anomaly 
fields,  from  which  can  be  determined  seasonally  averaged  geostrophic  current  velocity  fields.  The 
location  of  the  data  is,  however,  generally  temporally  and  spatially  random,  complicating  the 
averaging  procedure,  and  thus  requiring  the  utilization  of  special  methods. 

THE  METHOD 


The  most  direct  method  to  analyze  the  data  is  to  group  together  all  the  dynamic  depth  anomaly 
data  for  a  particular  season  and  then  fit  this  data  by  the  method  of  least  squares  with  a  two 
dimensional  polynomial  of  the  form: 

M  N 

dD  =  Z  2  *"  (1) 

m=o  n=o 

where 


do  =  dynamic  depth  anomaly  at  depth  D 

Amn  =  regression  coefficients 

d  =  latitude 

<p  -  longitude 

M,N  =  order  of  the  polynomials 

The  virtue  of  this  method  is  that  it  treats  each  data  point  for  exactly  what  it  is,  an  independent 
measure  of  the  dynamic  depth  field,  and  at  the  same  time  yields  an  analytic  expression  for  the 
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dynamic  depth  anomaly  field.  The  latitude  and  longitude  are  used  as  the  independent  variables 
because  this  is  the  most  convenient  form  in  which  to  obtain  the  data  from  N ODC  and  thus  no 
coordinate  transformations  are  necessary.  The  north-south  and  east-west  geostrophic  velocities 
can  now  be  obtained  by  differentiating  "  /ration  (1): 

Ufl  =  10  9(i  D  (2) 

f  R  sin  9  3  <p 

v .  = _ ll _  iia 

v  v  f  r  be 


where 

U$  =  north-south  velocity 

V  ^  =  east-west  velocity 

f  =  Coriolis  Parameter 

R  =  Radius  of  the  earth 


The  dynamic  depth  anomaly  field  frequently,  however,  displays  so  much  structure  that  to 
replicate  it  accurately  with  a  single  function  of  the  form  given  in  equation  (1)  would  require  a 
polynomial  of  prohibitively  high  order.  This  difficulty  can  be  overcome  through  the  use  of  two- 
dimensional  spline  fits  of  the  data.  In  this  case  the  area  of  interest  is  divided  into  rectangles 
and  the  data  in  each  rectangle  is  still  fit  by  the  method  of  least  squares  with  a  polynomial  of  the 
form  expressed  in  equation  (1),  but  with  the  constraints  that  the  different  polynomials  and  their 
first  and,  if  necessary,  higher  derivative  match  on  the  boundaries  of  the  rectangles.  The  order  of 
the  polynomial  is  selected  to  yield  whatever  goodness  of  fit  is  desired.  The  advantage  of  spline 
fitting  is  that  the  data  which  was  previously  fitted  by  one  polynomial  is  now  fitted  by  several 
polynomials,  allowing  much  more  structure  to  be  represented.  Equations  (2)  can  still  be  used  to 


evaluate  the  components  of  geostrophic  velocity. 
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THE  PH  OCR  AM 


To  put  this  method  into  use  a  spline  fit  program  called  SPLPAK  which  employs  third  order 
polynomials  was  obtained  form  the  National  Center  for  Atmospheric  Research  (N'CAR)  Software 
Support  Library.  The  SPLPAK  program  has  several  options,  the  most  important  of  which  is 
the  number  of  rectangles  into  which  the  area  to  be  fitted  is  to  be  divided.  Since  it  was  assumed 
that,  even  with  spline  fits,  the  area  to  be  analyzed  would  be  too  large  to  be  fitted  by  a  single 
spline  fit,  it  was  specified  for  the  program  to  have  simply  nine  rectangles,  as  shown  in  Figure  1, 
with  the  size  of  the  rectangles  being  left  as  a  variable  to  be  determined  bv  the  needs  of  the  parti¬ 
cular  area  under  study.  The  idea  behind  using  nine  rectangles  was  that  the  best  fit  to  the  data 
should  occur  in  the  center  rectangle  where  the  effects  of  the  outer  boundary  are  least  significant. 
The  total  area  of  interest  would  then  be  analyzed  by  chosing  an  appropriate  size  for  the  nine 
rectangles  according  to  the  amount  of  structure  present  and  then  moving  the  center  rectangle 
around  the  total  area  of  interest  with  some  continuity  being  supplied  from  one  position  of  the 
center  rectangle  to  an  adjacent  position  by  means  of  the  overlapping  effect  of  the  eight  sur¬ 
rounding  rectangles.  This  approach  would  allow  different  size  areas  to  be  analyzed  with  equal 
accuracy  and  also  allow  for  the  development  of  a  generalized  program  applicable  to  any  area. 
A  fuller  explanation  of  SPLPAK  and  its  other  options  can  be  found  in  Adams  et  al.  (1975). 

In  order  to  utilize  SPLPAK  a  program  called  XTRACT  was  written  to  read  a  standard  NODC 
data  tape  and  write  on  disk  the  appropriate  input  data  required  by  SPLPAK.  A  program  called 
DYF1T  was  then  written  which  would  read  the  data  from  the  desk,  call  SPLPAK  as  a  subroutine 
and  evaluate  the  components  of  the  geostrophic  velocity  at  any  specified  location.  Besides  the 
size  of  the  rectangles,  the  relative  depth  to  wiiich  the  currents  are  to  be  referenced  and  the  time 
period  over  which  the  currents  are  to  be  averaged  are  variables  which  must  be  specifier!  in  DYFIT. 


A  listing  and  fuller  explanation  of  the  programs  XTRACT  and  DYF1T  can  he  found  in  Appendices 
A  and  B. 

THE  SOUTHERN  CALIFORNIA  BIGHT 

A  region  from  32°  N  to  3a  N  and  from  1 1 T  ^  to  122°W  off  the  coast  of  Southern  California 
was  chosen  as  the  area  to  test  the  applicability  of  the  above  described  method.  There  were  several 
reasons  for  choosing  this  particular  area: 

o  There  have  been  many  hydrographic  surveys  in  the  area,  so  there  should  be  a  sufficient 
amount  of  data: 

o  Qualitatively  the  currents  arc  well  understood  showing  a  definite  seasonal  dependence: 

o  The  water  is  deep  so  that  the  geostrophic  currents  should  be  meaningful. 

A  standard  hydrographic  data  tape  for  the  above  area  w’as  obtained  from  NODC.  After  initial 
preparation  of  the  data  with  the  XTRACT  program  the  DYFIT  program  was  run  with  rectangles 
equal  to  1°  of  latitude  by  1°  of  longitude.  Also  as  most  of  the  hydrograpluc  data  extended  only  to 
500  meters,  500  meters  was  chosen  as  the  reference  level,  thus  enabling  the  utilization  of  most 
of  the  data.  While  500  meters  may  appear  to  be  somewhat  shallow  for  a  reference  level,  cross¬ 
checking  wherever  possible  showed  no  significant  differences  in  surface  velocities  when  1000 
meters  was  used  as  a  reference  level.  The  spline  fit  was  evaluated  at  the  center  of  rectangles  1/4°  x 
1/4°  1/2°  x  1/2°  and  1°  x  1°  on  their  sides.  As  it  turned  out  the  current  structure  was  such  that 
there  were  no  areas  that  required  rectangles  smaller  than  1/2°  x  1/2°  in  order  to  be  able  to  cha¬ 
racterize  the  currents.  The  data  was  of  sufficient  volume  to  be  grouped  by  months  for  analysis. 

The  results  of  the  analysis  are  depicted  in  the  twelve  chartlets  in  Figures  2-13  with  the  current 
directions  given  in  degrees  true  and  the  speed  in  knots.  The  currents  are  evaluated  at  the  center 
of  the  rectangle.  The  analysis  indicates  that  the  average  geostrophic  velocities  in  the  Southern 
California  Bight  are  relatively  small,  with  maximum  speeds  of  onlv  approximately  0.2  knot. 
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Unfortunately  an  in-depth  critical  quantitative  comparison  of  observed  and  predicted  currents  is 
not  possible  due  to  a  lack  of  detailed  current  measurements  in  the  area.  The  calculated  geos- 
trophic  velocities  do,  however,  conform  with  the  accepted  general  circulation  scheme  for  the  area 
as  presented  for  example  in  Pirie  and  Stellar  (1977).  This  scheme  presents  the  Southern  California 
Bight  as  being  dominated  by  a  cyclonic  gyre,  with  the  southeastwardlv  flowing  California  Current 
to  its  west,  the  northwestwardly  flowing  Davidson  Current  to  its  east.  The  speeds  in  the  gyre  are 
known  to  be  seasonally  dependent.  On  the  west  the  current  speeds  increase  during  the  spring  and 
early  summer  as  the  California  Current  tends  to  move  eastward  during  this  period.  The  Davidson 
Current  is  an  undercurrent  during  this  time  period  with  little  surface  manifestation.  During  the 
fall  and  winter  the  Davidson  Current  surfaces  and  the  California  moves  westward  with  a  subse¬ 
quent  increase  in  the  current  speed  on  the  eastern  side  of  the  gyre  and  a  decrease  of  those  on  the 
west.  While  it  is  realized  that  rigorous  comparisons  are  not  possible  without  precise  measurements, 
the  above  described  scheme  is  well  supported  qualitatively  in  the  tw  elve  chartlets  in  Figures  2-13. 
CONCLUSIONS 

A  method  has  been  developed  which  uses  two-dimensional  least-squares  spline  fits  of  spatially 
and  temporally  random  dynamic  depth  anomaly  data  to  determine  average  geostrophie  currents. 
The  method  has  successfully  been  applied  to  quantify  the  currents  in  the  Southern  California 
Bight.  It  can  be  assumed  that  it  is  possible  to  apply  this  method  to  other  areas  and  thus  begin  to 
utilize  the  vast  amounts  of  data  in  the  NODC  inventory . 
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APPENDIX  A 


I.  Propram  Name:  XTRACT 

II.  System  Name.  DYFIT 

III.  Purpose:  XTRACT  is  the  first  step  in  a  system  of  proprams  that  build  a  data  file  for 

the  DYFIT  propram.  XTRACT  reads  NO.AA  STD  Data  tape  and  extracts  records  with  valid 
dynamic  depth  anomaly  data.  All  other  records  are  skipped.  The  data  is  processed  (see  section 
VII)  and  written  out  onto  a  disk  file  for  subsequent  processing  by  the  IBM  SORT  routine. 

IV.  History:  XTRACT  was  originally  written  in  FORTRAN  by  David  A.  Portyrata  in  October 
1978.  It  was  rewTitten  in  COBOL  in  July  1979.  The  change  to  COBOL  allowed  the  use  of  built- 
in  functions  to  handle  the  overpunched  sign  on  the  latitude  and  longitude  fields,  and  facilitated 
use  of  packed  decimal  output. 

V.  Machine:  IBM  360/65 

VI.  Language:  IBM  COBOL  VERSION  4 

VII.  Description  of  Processing: 

a.  Scanning  for  Valid  Records 

Many  different  types  of  data  records  are  contained  on  a  STD  tape.  Only  record 
type  6  contains  a  dynamic  depth  anomaly  field.  All  other  records  are  skipped.  The  following 
conditions  are  necessary  for  a  record  to  be  processed: 

(1)  A  '6'  in  the  80th  byte  (the  record  type  field). 

(2)  Valid  numeric  data  in  bytes  60-63  (the  dynamic  depth  anomaly  field). 

(3)  Valid  numeric  data  in  bytes  7-8  (the  minutes  of  latitude  field). 

(4)  Valid  numeric  data  in  bytes  13-14  (the  minutes  of  longitude  field). 

b.  Converting  Latitude  and  Longitude  Records 

Latitude  and  longitude  records  are  converted  from  degrees/minutes/tenths 
of  minutes  to  decimal  degrees.  The  latitudes  are  subtracted  from  90  so  that  they  range  from 
0  to  180  instead  of  -90  to  90.  180  is  added  to  the  longitudes  so  that  they  range  from  0  to  360 
instead  of  -180  to  180.  The  above  operations  are  performed  so  that  the  output  from  XTRACT 
need  not  contain  signed  numeric  fields. 

c.  Computing  the  Julian  Day 

The  Julian  date  is  calculated  using  a  formula  published  in  the  Texas  Instru¬ 
ments  Master  Library'  for  TI  58/59.  The  factor  for  the  cast  day  is  calculated,  and  the  factor  for 
January  1st  of  the  year  the  cast  was  taken  is  subtracted  from  it,  giving  the  Julian  date.  The 
formula  is:  A , 


For  January  and  February: 


FACTOR  =  365  (YYYY)  +  DD  +  31  (MM  -  1)  +  INT(YYYY  -  1)1  A) 

■  INT(3/4[INT]  ((YYYY  -  1)/100)  +  1]) 

For  March  through  December: 

FACTOR  =  365  (YYYY)  +  DD  +  31(MM  -  1)  -  INT(4MM  +  23) 

+  INT  (YYYY/4)  ■  INT(3/4[INT  (YYYY/100)  +  1]) 

In  these  formulas  INT  indicates  using  only  the  integer  portion  of  the  expression. 

The  number  of  days  between  two  dates  is  found  by  calculating  the  FACTOR  for  each 
date,  then  finding  the  difference  between  the  FACTORS. 

VIII.  File  Description: 


a.  Input  File 

Input  to  XTRACT  is  a  NOAA  tape  of  STD  data.  Specific  information  as  to 
tape  characteristics  will  be  on  the  tape  itself.  The  following  characteristics  are  assumed: 


TRACKS  9 

DENSITY  1600  Bits/Inch 

LOGICAL  RECORD  LENGTH  80  Bytes 

BLOCKSIZE  3200  Bytes 

RECORD  FORMAT  FB 

UNLABELED 


b.  Output  File 

Output  from  XTRACT  is  a  disk  file  in  packed  decimal  form.  It  is  to  read 
into  the  IBM  SORT  utility.  The  following  characteristics  are  used: 

LOGICAL  RECORD  LENGTH  18  Bytes 

BLOCKSIZE  6372  Bytes 

c.  Structures  Used 

The  structures  used  for  input  and  output  of  data  in  the  EXTRACT  program 
are  shown  in  Tables  A1  and  A2. 


TABLE  A1 


STD  DATA  FILE  STRUCTURE  (XTR1N) 
(Logical  Record  Length  =  80  Bytes) 


POSITION 

DESCRIPTION 

01-04 

Filler 

05-09 

Latitude  -  positions  5  and  6  degrees,  positions  7  and  8 
minutes,  position  9  tenths  of  minutes,  eleven  overpunch 
in  position  8  for  south  latitude. 

10-15 

Longitude  -  positions  10-12  degrees,  positions  13  and  14 
minutes,  position  15  tenths  of  minutes,  eleven  overpunch  in 
position  14  for  east  longitude 

16-18 

Filler 

19-24 

Date  -  positions  19  and  20  years,  positions  21  and  22  month 
and  positions  23  and  24  day. 

25-27 

Filler 

28-31 

Standard  Depth 

32-59 

Filler 

60-63 

Dynamic  depth  anomaly 

64-79 

Filler 

80 


Record  type  (6) 


TABLE  A2 


STANDARD  DEPTH  FILE  STRUCTURE  (XTROUT) 
(Logical  Record  Length  =  18  Bytes) 


POSITION 

DESCRIPTION 

01-03 

Standard  depth 

04&05 

Julian  day 

05&07 

Year 

08-11 

Latitude-decimal  degrees 

12-15 

Longitude-decimal  degrees 

16-18 

Dynamic  depth  anomaly 

A4 


IDENTIFICATION  DIVISION. 

PROGRAM- IO.  XTRACT. 

AUTHOR.  BOR  STAR*.  CO-OP.  WEST  VIRGINIA  INSTITUTE  OF  TECHNOLOGY. 
INSTALLATION.  US  COAST  GUARO  OCEANOGRAPHIC  UNIT.  ATD/DP9. 
DATE-WRITTEN.  JUL  2.  1980. 

DATE-COMPILED. 

REMARKS. 

THIS  PROGRAM  SCANS  NODC  STD  TAPE  RECORDS  FOR  STANDARD 
DEPTH  DYNAMIC  OEPTH  DATA.  WHEN  A  STANDARD  DEPTH  IS  FOUND. 

THE  JULIAN  DAY  FOR  THAT  RECORD  IS  CALCULATED.  THEM  THE 
STANDARD  DEPTH,  JULIAN  DAY,  YEAR,  LATITUDE,  LONGITUDE.  AND 
DYNAMIC  DEPTH  ARE  WRITTEN  OUT  ONTO  A  FILE  IN  PACKED  DECIMAL 
FOPMAT.  THIS  DATA  WILL,  AFTER  SUBSEQUENT  PROCESSING  RY  CREATE 
AMD  DUPFIX.  RECODE  INPUT  TO  THE  D YE  I T  PROGRAM. 

THIS  RROGRAM  IS  A  REWRITE  OF  THE  ORIGINAL  FORTRAN  XTRACT 
PROGRAM,  WRITTEN  OCT  2,1979  RY  DAVID  PORTYRATA.  IT  WAS 
REWRITTEN  IN  COROL  TO  SIMPLIFY  THE  HANDLING  OF  THE  OVERPUNCH 
MINUS  SIGN  IN  THE  LATITUDE  AND  LONGITUDE  FIELDS. 

*  NOTE  *...Th£  JULIAN  DAY  IS  CALCULATED  USING  AN  ALGORITHM 
FROM  "T . I .  PROGRAMMABLE  5B/S9  MASTER  LIBRARY",  PAGE  76. 

ENVIRONMENT  DIVISION. 

CONFIGURATION  SECTION. 

SOURCE-COMPUTER.  IBM-360. 

OBJECT-COMPUTER.  IBM-3H0. 

INPUT -OUTPUT  SECTION. 

FILE-CONTROL. 

SELECT  STD-DATA-EILE  ASSIGN  TO  UT-S-XTRIN. 

SELECT  STANDARD-OEPTH-EILE  ASSIGN  TO  UT-S-XTROUT. 

DATA  DIVISION. 

FILE  SECTION. 

FD  STO-DATA-FILE  LABEL  RECORDS  OMITTED 
RECORD  CONTAINS  80  CHARACTERS 
BLOCK  CONTAINS  0  RECORDS. 

01  STD-DATA-REC. 


03 

FILLER 

PIC 

X  (  A)  . 

03 

LATITUDE. 

05 

LAT-OEGREES 

PIC 

9(2) 

05 

lat-minutes 

PIC 

S9  ( 2  > 

05 

LAT-1 OTH-MIN 

PIC 

9. 

03 

LONGITUDE. 

05 

LON-OEGREES 

PIC 

9(3) 

05 

lon-minutes 

PIC 

S9  (?) 

05 

LON-1 OTH-MIN 

PIC 

9. 

03 

FILLER 

PIC 

X  ( 3 )  . 

03 

DATE 

-DATA. 

05 

DATE-YEAR 

PIC 

9(2). 

05 

DATE-MONTH 

PIC 

9(2)  . 

05 

OATE-OAY 

PIC 

9(2)  . 

03 

filler 

Pic 

X  (3)  . 

03 

STANDARD-DEPTH 

PIC 

9(A)  . 

03  FILLER  PTC  X(?R). 

03  DYNAMIC-DEPTH  PIC  S9V999. 

03  FILLE°  PIC  X(lft). 

03  STANDARD-DEPTH-FLAG  PIC  9. 

FD  STANQA9D-OEPTH-FILE  LAPEL  RECORDS  STANDARD 
PECOPO  CONTAINS  IP  CHAPACTFPS 
BLOCK  CONTAINS  354  RECORDS. 

01  STANDARD-DEPTH-REC  USAGE  COMP-3. 


03 

ST-STANDARD-DFPTH 

PIC 

9  (S)  . 

03 

ST-JULI AN-DAY 

PIC 

9(3). 

03 

ST-YEAR 

PIC 

9(2)  . 

03 

st-latitude 

PIC 

999V999. 

03 

ST-LONGITUDE 

PIC 

999V999. 

03 

ST-DYNAMIC-DEPTH 

PIC 

S99V999. 

WORKING-STOWAGE  SECTION. 

77  EOF-FLAG 
77  MINUTES 
77  REC-COUNT 
77  READ-COUNT 
77  WS-YEAR 
77  WS-LATITUDE 
77  WS-LONGITUDE 
77  REAL-1 
77  REAL-2 
77  REAL-3 
77  INTEGER-1 
77  INTEGER-2 
77  INTEGER-3 
77  JAN- 1 -FACTOR 
77  FACTOR 

PROCEDURE  DIVISION. 

*  THIS  MODULE  IS  THE  SUPERVISOR.  HOUSEKEEPING  IS  DONE  HERE. 
OPEN  INPUT  STO-DATA-FILE 

OUTPUT  STANDARD-DEPTH-FILE. 

PERFORM  WEAO-REC 

PERFORM  MAIN-PROCESS  THRU  MAIN-PROCESS-EXIT 
UNTIL  EOF-FLAG  =  1. 

DISPLAY  ♦  END  OF  PROCESSING  BY  XTRACT  *.  READ-COUNT  * 

•  RECORDS  READ  «♦  REC-COUNT.  •  RECORDS  WRITTEN.*. 
CLOSE  STD-DATA-FILE 

standard-depth-file. 

STOP  RUN. 


MAIN-PROCESS. 

*  THIS  MODULE  MAKES  SURE  THE  DATA  IS  A  STANOARO  DEPTH. 

*  AND  CONVERTS  THE  LATITUOE  AND  LONGITUDE  TO  THE  PROPER  FORM. 
ADO  1  TO  READ-COUNT. 

EOIT-CHECK-THE-OATA. 

IF  LAT-OEGREFS  IS  NOT  NUMERIC 


PIC  9  USAGE  COMP  SYNC. 

PIC  99  USAGE  CO*P  SYNC. 

PIC  9(B)  VALUE  0  USAGE  COMP  SYNC. 

PIC  9(B)  VALUE  0  USAGE  COMP  SYNC. 

PIC  9(B)  USAGE  COMP  SYNC. 

PIC  S999V999. 

PIC  S999V999. 

USAGE  COMP-1  SYNC. 

USAGE  COMP-1  SYNC. 

USAGE  COMP-1  SYNC. 

PIC  9(B)  USAGE  COMP  SYNC. 

PIC  9(B)  USAGE  COMP  SYNC. 

PIC  9(B)  USAGE  COMP  SYNC. 

PIC  59 ( P )  USAGE  COMP  SYNC. 

PIC  S9(B>  USAGE  COMP  SYNC. 


OR  LON-DEGREES  IS  NOT  NUMERIC 
OR  DATE-DATA  IS  NOT  numeric 
OR  STANDARD-DEPTH  is  NOT  NUMERIC 
op  STANDAoD-OEPTH-ELAG  TS  NOT  NUMERIC 
Then  PPREORM  REaD-REC 

GO  TO  MAIN-RROCESS-EX IT  . 

CHECK  FOR  A  STANDARD  DEPTH.  IF  NOT.  READ  NEXT  REC.  AND  EXIT. 
IE  STANDARQ-OEPTH-ELAG  NOT  =  6 

OR  DYNAMIC-DEPTH  IS  NOT  NUMERIC 
THEN  PEPFORm  READ-RFC 

GO  TO  MAIN-PROCESS-EXIT . 

convert  degrcfs-minutes-tenth-minutes  OF  latitude  TO  DECIMAL 
DEGREES.  (NOTE.. .MINUTES  IS  AN  UNSIGNED  FIELD.) 

IF  LAT-MINUTFS  IS  NOT  numeric 
THEN  PERFORM  READ-REC 

GO  TO  maIN-PROCESS-EXIT 
ELSE  MOVE  LAT-MINUTFS  TO  MINUTES. 

CONVERT  BLANKS  TO  ZERO. 

IF  LAT-IOTH-MIN  IS  NOT  NUMERIC 
Then  MOVE  0  TO  LAT-1 OTH-MIN. 

COMPUTE  WS-LATITUDE  =  LAT-OEGREES  ♦  (MINUTES  /  GO) 

♦  (LAT-1 OTH-MIN  /  600)  . 

IE  LAT-MINUTFS  IS  NEGATIVE 

THEN  MULTIPLY  WS-LATITUDE  BY  -1  GIVING  WS-LATITUDE. 

90  IS  ADDED  TO  ALL  LATITUDES  SO  THEY  WILL  BE  POSITIVE 
AND  RAN6E  FROM  o  TO  ISO. 

ADD  90.  WS-LATITUDE  GIVING  ST-LATITUDE. 

CONVERT  DEGREES-MINUTES- TENTHS  OF  LONGITUDE  TO  DECIMAL  DEGREE 
IE  LON-MINUTFS  IS  NOT  NUMERIC 
THEN  PERFORM  REAO-REC 

GO  TO  maIN-PROCESS-EXIT 
ELSE  MOVE  LON-MINUTES  TO  MINUTES. 

CONVERT  BLANKS  TO  7ER0. 

IF  LOM-IOTH-MIN  TS  NOT  NUMERIC 
THEN  MOVE  0  TO  LON-1 OTH-MIN. 

COMPUTE  WS-LDNGITUDE  =  LON-DEGPEES  ♦  (MINUTES  /  GO) 

♦  (LON-10TH-MIN  /  600). 

IF  LON-MINUTFS  IS  NOT  NEGATIVE 

THEN  MULTIPLY  WS-LOUGI TUOE  BY  -1  GIVING  WS-LONG I TUDE . 

I H 0  IS  AOOEO  TO  ALL  LONGITUDES  SO  THAT  THEY  WILL 
BE  POSITIVE  AND  PANGE  FROM  0  TO  160. 

ADD  1  B  0 .  wS-LONGITUDE  GIVING  ST-LONG  I TUDE . 

COMPUTE-JULI AN-OAY. 

FIRST  WE  FIND  THE  JULIAN  DATE  FACTOR  FOR  JANUARY  1ST  OF 
THE  YEAR  THE  CAST  WAS  TAKEN. 

ADD  1900.  DATE-YEAR  GIVING  WS-YEAR. 

COMPUTE  INTEGER-1  =  (WS-YEAR  -  ])  /A. 

COMPUTE  INTEGER-?  =  < (WS-YEAR  -  l)  /  100)  *  1 
COMPUTE  REAL-?  =  INTEGER-?  *  0.76. 

MOVE  REAL-?  TO  INTEGER-?. 

COMPUTE  JAN-1-FaCTOo  =  (365  *  WS-YEAR)  ♦  INTEGER-1  ♦  1  - 


INTEGER-? 


*  NOW  w£  FIND  THE  JULIAN  DATE  FACTOR  FOR  THE  DAY  THE  CAST 

*  WAS  TAKEN. 

IF  DATE-MONTH  <  3 

then  perform  jan-fer-factor 
ELSE  PERFORM  MaR-OEC-FACTOR. 

*  NOW  WE  FIND  THE  DIFFERENCE  IN  THE  FACTORS*  THE  JULIAN  DAY 
COMPUTE  ST-JULI AN-DAY  =  FACTOR  -  JAN-l-FACTOR  ♦  1. 

*  MOVE  The  kEST  OF  THE  FIELDS  TO  THE  OUTPUT  BUFFER. 

MOVE  DATE-YEAR  TO  ST-YEAR. 

MOVE  STANDARD-DEPTH  TO  ST-ST ANDARD-DEPTH. 

MOVE  DYNAMIC-DEPTH  TO  ST-QVNAMIC-DEPTH, 

*  WRITE  OUT  THF  RECORD  A  READ  IN  DATA  FOR  THE  NEXT  LOOP. 

WRITE  STANPARO-DEPTH-HEC. 

ADD  1  TO  REC-COUNT. 

PERFORM  WEAD-REC. 

MAIN-RROCESS-EXIT. 

EXIT, 

JAN-FEB-FACTOR. 

*  this  module  computes  Julian  date  factors  for  January  and 

*  FEBRUARY.  ALGORITHM  FROM  TEXAS  INSTRUMENTS  MASTER  LIBRARY. 

«  NOTE... INTEGER-1  &,  2  WERE  ALREADY  CALCULATED  WHEN 

*  JAN-l-FACTOR  WAS  calculated. 

COMPUTE  FACTOR  =  t 3G5  *  wS-YE AR )  ♦  (31  *  (DATE-MONTH  -  p 
)  ♦  DATE-DAY  ♦  INTEGER-1  -  INTEGER-2. 

“AR-OEC-FACTPR. 

*  THIS  MODULE  COMPUTES  JULIAN  DATE  FACTORS  FOR  MARCH  THRU 

*  DECEMBER.  ALGORITHM  FROM  TEXAS  INSTRUMENTS  MASTER  LIBRARY. 

COMPUTE  HEAL-1  =  (0.4  *  DATF-MONTH)  ♦  2.3. 

MOVE  RFaL-1  TO  INTEGER-1, 

CUMRuTF  INTEGER-?  =  WS-YEAR  /  4. 

COMPUTE  INTEGER-3  =  (wS-YEAR  /  100)  ♦  1. 

COMPUTE  REAL-3  =  INTEGER-3  *  0.75. 

MOVE  REaL-3  TO  INTEGER-3. 

COMPUTE  FACTOR  a  ( 3*5  *  WS-YEAR)  ♦  (31  *  (DATE-MONTH  -  1) 

)  ♦  OATE-OAY  -  INTEGER-1  ♦  INTEGER-2  -  I NTEGER-3 . 

READ-PEC. 

*  THIS  MODULE  READS  RECORDS  FROM  THE  STD  DATA  TAPE. 

READ  STD-DATA-FILE  AT  END  MOVE  1  TO  EOF -FLAG. 


APPENDIX  B 

I.  Program  Name:  DYFIT 

II.  System  Name:  DYFIT 

III.  Purpose-.  The  DYFIT  program  calculates  the  average  geostrophic  velocity  field  from 
randomly  located  dynamic  depth  anomaly  data. 

IV.  History:  The  DYFIT  program  was  written  by  K.  A.  Mooney  for  the  CDC  3300  computer 
in  the  summer  of  1978.  Due  to  lack  of  memory  and  slow  processing  speed  of  the  CDC,  this 
version  was  impractical  for  production  use.  The  IBM  360  version  was  written  by  K.  A.  Mooney 
in  the  spring  of  1979.  This  version  was  used  to  study  currents  in  the  Southern  California,  Gulf 
of  Mexico,  and  Hawaii  regions.  This  version  searched  a  large  sequential  file  (created  by  XTRACT 
or  XTEMP)  for  suitable  dynamic  depth  anomaly  data.  This  program  wras  highly  I/O  bound 
because  it  only  used  about  0.5%  of  the  records  that  it  has  to  read.  The  current  version  of  DYFIT 
was  written  by  Bob  Stark  in  the  summer  of  1980.  This  version  utilizes  a  COBOL  subroutine 
which  reads  small  sections  of  an  indexed  sequential  file.  This  reduced  the  I/O  time  required  to 
run  the  program.  The  XTRACT  and  XTEMP  programs  were  rewritten  and  the  CREATE  and 
DUPFIX  programs  written  to  build  indexed  sequential  files  for  this  version. 

V.  Machine:  IBM  360/65 

VI.  Language: 

a.  Languages  Used 

1.  IBM  ‘H’  level  FORTRAN  (Optimized) 

2.  IBM  version  4  ANSI  COBOL  (Optimized) 

b.  Details  of  FORTRAN-COBOL  linkage 
1.  Compilation 

The  COBOL  modules  (DYFIT  and  SETKEY /READER)  are  compiled  in  separate 
runs  of  the  COBOL  compiler.  The  object  module  produced  by  the  first  run  (ddname  =  SYSLIN) 
is  passed  on  to  the  second  run.  All  of  the  FORTRAN  subroutines  are  compiled  in  a  single  run  of 
the  FORTRAN  compiler  and  placed  in  the  same  data  set  as  the  COBOL  object  modules  (SYSLIN). 
The  linkage  editor  is  then  run  to  link  all  of  the  object  modules  into  a  load  module  and  store  it  on 
disk.  The  details  of  this  process  may  be  seen  in  the  cataloged  procedure  K198PDFT,  which  is 
used  to  put  the  DYFIT  program  to  disk. 

Transfer  of  control  from  main  program  to  subroutine  and  back  is  accomplished 
by  the  DYFIT  program,  which  calls  a  ILBOSTPO  (a  library  subroutine)  to  initialize  the  sub¬ 
routine  communications  save  area.  This  assures  that  all  routines  called  by  DYFIT,  or  called  by 


routines  called  by  DYFIT,  will  be  treated  as  subroutines.  This  is  especially  important  for  the 
SETKEY/READER  routine,  as  EXIT  PROGRAM  statements  in  the  subroutine  communications 
save  area  in  order  for  control  to  return  to  the  calling  routine. 

3.  Correspondence  of  Identifiers 

Two  types  of  identifiers  were  used  as  parameters  in  the  COBOL  routines  to  be 
called  by  FORTRAN  routines:  Fullword  real  (floating  point)  and  halfword  integer  (binary). 
These  variables  must  be  aligned  on  fullword  boundaries,  which  is  automatic  in  FORTRAN.  In 
COBOL  the  keyword  SYNCRONIZED  (abbreviated  SYNC)  is  used  to  align  the  identifiers  on  a 
fullword  boundary. 

The  fullword  real  identifiers  are  declared  as  R  F.AL*4  in  the  FORTRAN  programs, 
and  defined  as  USAGE  COMP-1  SYNC  in  the  COBOL  routines.  These  two  definitions  are  fully 
compatable. 

The  halfword  integer  identifiers  are  declared  as  INTEGER*  2  in  the  FORTRAN 
programs,  and  defined  as  PIC  S9999  USAGE  COMP  SYNC  in  the  COBOL  routines. 

VII.  Description  of  Processing:  DYFIT  calls  DYFORT  which  reads  a  parameter  card  specifying 
a  latitude,  longitude,  and  time  range.  This  card  also  indicates  from  what  standard  depth  the  data 
is  to  come.  SF.TKEY  is  called  to  position  the  indexed  sequential  file  at  the  first  record  of  the 
user  specified  standard  depth  that  falls  within  the  time  range  chosen.  The  file  is  searched  from 
this  point  for  dynamic  depth  anomaly  data  that  falls  within  the  specified  latitude  and  longitude. 
The  search  stops  when  the  given  time  range  is  exceeded,  the  given  standard  depth  is  exceeded, 
or  the  physical  end  of  file  is  reached.  SPLCW  is  then  called  to  fit  a  two  dimensional  cubic  spline 
to  the  dynamic  depth  anomaly  data.  A  description  of  SPLCW  can  be  found  in  the  chapter  on 
SPLPAK  in  volume  1  of  the  NCAR  SOFTWARE  SUPPORT  LIBRARY,  NCAR  technical  note 
IA-105,  March,  1975.  Finally,  DYFORT  evaluates  the  first  derivatives  of  the  spline  fit  at  user 
specified  locations  by  calling  the  function  SPLDE.  A  description  of  SPLDE  can  also  be  found  in 
the  above  reference.  The  components  of  the  geostrophic  velocities  can  then  be  determined  from 
the  formulae: 


where 

Ug  =  North  -  South  component  of  geostrophic  velocity" 
V  ^  =  East  -  west  component  of  geostrophic  velocity 


f  =  Coriolis  parameter 

R  =  Radius  of  earth 

dp  =  Dynamic  depth  anomaly  field  at  depth  D 
0  =  Longitude 

0  =  Latitude 

The  printout  consists  of  the  standard  depth  of  the  data  and  the  center  and  ranges  for  the 
latitude,  longitude,  and  Julian  day.  At  the  end  of  the  data  list  (appears  if  data  print  option  was 
selected),  the  latitude,  longitude,  and  Julian  day  at  which  the  plane  fit  was  evaluated  along  with 
the  components  of  geostrophic  velocity  in  knots  and  the  number  of  points  which  were  used  in 
the  fit  are  printed. 

VIII.  File  Description: 

a.  Disk  File 

1.  Physical  Characteristics 

(a)  Indexed  sequential  organization 

(b)  Sequential  access  mode 

(c)  One  cylinder  used  for  cylinder  index 

(d)  RECFM  =  FB  (fixed  longth  blocks) 

(e)  LRECL  =  18  bytes  (logical  record  length) 

(f)  KEYLENGTH  =  15  bvtes 

(g)  BLKSIZE  =  6372  bytes 

(h)  Blocking  factor  354  records/block 

(i)  All  fields  in  packed  decimal  form. 

2.  Usage 

The  dynamic-height-file  (ISAMFILE)  contains  the  dynamic  depth  anomaly  data 
for  DYFIT.  It  is  created  by  the  CREATE  program  -  refer  to  it  for  sort  order.  In  DYFIT,  a  call 
to  SETKEY  opens  the  file  and  positions  it  at  the  first  record  matching  the  starting  -  standard  - 
depth  (DPDP)  and  the  first  -  Julian  -  day  (TMIN)  parameters.  If  there  is  no  data  in  the  file  at 
this  position,  SETKEY  returns  a  1  in  key  -  error  (KEYERR).  READER  is  then  called  to  read 
records  sequentially  from  that  point  on  until  a  record  is  read  whose  Julian  day  is  greater  than 
last  -  Julian  -  day  (TMAX),  whose  standard  -  depth  is  greater  than  starting  -  STD  -  depth  (DPDP), 
or  the  end  of  file  is  found.  The  file  is  then  closed  and  reader  returns  a  1  in  EOF. 


b.  Card  File 


1.  Physical  Characteristics 

(a)  Sequential  organization 

(b)  LRECL  =  80  bytes  (logical  record  length) 

(c)  Unblocked 

c.  Structures  Used 

1.  The  structures  used  for  the  input  files  are  shown  in  Tables  Bl,  B2  and  B3. 


B4 


TABLE  B1 


DYNAMIC  HEIGHT  FILE  (1SAMFILE) 


POSITION 

DESCRIPTION 

01-03 

Standard  depth 

04&05 

Julian  day 

06&07 

Year 

08-11 

Latitude 

12-15 

Longitude 

16-18 

Dynamic  depth 
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TABLE  B2 


PARAMETER  CARD  (SYSIN) 


POSITION 

DESCRIPTION 

01-05 

Standard  depth 

06-08 

Filler 

09-14 

Center  latitude 

15 

Filler 

16-22 

Center  longitude 

23 

Filler 

24-27 

Center  Julian  day 

28 

Filler 

29-33 

Latitude  range 

34 

Filler 

35-40 

Longitude  range 

41 

Filler 

4245 

Julian  day  range 

46 

Filler 

47 

Print 

48-80 

Filler 
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TABLE  B3 


POSITION 

01-06 

"7 

( 

08-14 

15-80 


DATA  CARD  (SYSIN) 
DESCRIPTION 
Latitude 
Filler 
Longitude 
Filler 
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IDENTIFICATION  DIVISION. 

PMOGMaM-IU.  0  V  E  I  T  . 

AUTHOR.  BOR  STARK,  CO-OP.  WEST  VIRGINIA  INSTITUTE  OF  TECHNOLOGY. 
INSTALLATION.  US  COAST  GUARD  OCEANOGRAPHIC  UNIT.  ATD/DPH. 

0  A  T  E  -  w  R  I  T T F N .  J'JN  ?R,19M0. 

DATE-COMPILED. 

REMARKS.  This  PROGRAM  IS  The  maim  PROGRAM  OF  The  DyFIT  SYSTEM. 

IT  HAS  P  STATEMENTS  -  A  CALL  TO  OYFQrT  (Th£  OLD  MAIN  POUTINE) 
AND  A  STOP  STATEMENT. 

The  first  COROL  Program  in  a  MIxFfl  LANGUAGE  ENVIRONMENT 
CALLS  TLROSTPO  (A  LlRRARY  SUBROUTINE)  WhICh  INITIALIZES  ThE 
SURROUTINF  COMMUNICATIONS  AREA  AND  SAVES  POINTER  REGISTERS 
NECESSARY  FOR  RETURN  TO  THE  MAIN  PROGRAM.  THEREFORE  ThE  FIRST 
PROGRAM  MUST  RF.  A  COROL  PROGRAM,  AS  CONTROL  WOULD  NOT  RETURN 
TO  A  MAIN  EO°TRAN  PROGRAM  AFTER  IT  HAD  CALLED  A  COROL  SUB¬ 
ROUTINE.  FOR  MORE  INFO.  REFER  TO  CALLING  AND  CALLED  PROGRAMS. 
IRM  OS  ANSI  COROL  PROGRAMMERS  GUIDE. 

ENVIRONMENT  DIVISION. 

DATA  DIVISION. 

PROCEDURE  DIVISION. 

CALL  • D Y F 0 R T  •  . 

STOP  RUN. 


IDENTIFICATION  DIVISION. 

PROGRAM- ID.  SETKFY. 

AUTHOR.  CO-OP  BOR  STARK  USCG  OCEANOGRAPHIC  UNIT. 

INSTALLATION.  US  COAST  GUARD  OCEANOGRAPHIC  UNIT.  ATD/OPR. 
DATE-wRITTFN.  *>/  1  1  /  1  PRO  . 

DATE-COMPILED. 

REMARKS.  this  PROGRAM  READS  AN  INDEXED  SEQUENTIAL  FILE 

CONTAINING  LATITUDE.  LONGITUDE.  YEAR.  JULIAN-DAY.  STANDARD 
DEPTH,  AND  DYNAMIC  HEIGHT  RECORDS;  THE  FILF  IS  SCWTED  RY 
STANDARD  DEPTH.  JULIAN  DAY,  LATITUDE.  AND  LONGITUDE ;  SKIP 
SEQUENTIAL  LOGIC  IS  USED  TO  POSITION  THE  FILE  AT  THE  PROPER 
START  point; 

TwO  ENTRY  POINTS  ARE  USED: 

1)  SETKFY  -  POSITIONS  THE  FILE  AT  ThE  FIRST  DESIRED 

RECORD  BASED  ON  THE  STD. DEPTH  AND  JULIAN  DAY 
?_)  READER  -  READS  THE  RECORDS  AND  RE  TU°NS  AN  END  OF 
FILE  FLAG*.  CLOSES  ThE  FILE  -HEN  E0F  =  1. 
ENVIRONMENT  DIVISION. 

CONFIGURATION  SECTION. 

SOURCE-COMPUTER.  lRM-3bO. 

OBJECT-COMPUTER.  IRm-3hq. 

INPUT-OUTPUT  SECTION. 

FILE-CONTROL. 

SELECT  OYNAMIC-HE IGHT-F ILE  ASSIGN  TO  0 A- I - I  SAMF I lE 
ACCESS  MODE  IS  SEQUENTIAL 
RECORD  KEY  IS  PfCORO-KEY. 


DATA  DIVISION 


file  section. 

FD  DYNAMIC-HEIGHT-FILE  LABEL  PECORDS  standard 
RECORD  CONTAINS  1*  CHARACTERS 
BLOCK  CONTAINS  354  RECORDS. 

01  DYNAMIC-HElGur-HFC. 

03  RECORD-xEY. 


OS  STANDARD-DEPTH  PIC 

9(5) 

USAGE 

COMP-3 

SYNC. 

OS  julian-day 

PIC 

RIG) 

USAGE 

COMP-3 

SYNC. 

OS  YEAR 

PIC 

9  (?) 

IJSAGF 

COMP-3 

• 

05  LATITUDE 

PIC 

999V99R 

USAGE 

COMP-3 

SYNC. 

OS  LONGITUDE 

PIC 

999V999 

USAGE 

COMP-3 

SYNC. 

03 

OYNAMIC-HE IGHT 

PIC 

S99V999 

USAGE 

COMP-3 

• 

WORKING- 

STORAGE  SECTION. 

01 

FOR- 

FLAG 

Pic 

S9999 

USAGF 

COMP 

SYNC. 

01 

GFNERIC-KEY. 

03 

NO'-'-STD-DERTh 

PIC  9 (S) 

USAGE 

COMP-3 

SYNC. 

03 

NQm-F IRST-JUL I  AN 

-DAY 

PIC  9(3) 

USAGE 

C  0  M  p  -  3 

SYNC. 

LINKAGE 

SECTION. 

& 

these  are  used  when 

setkey 

IS  CALLED. 

01 

F IRST-JUL I AN-OAY 

USAGE 

COMP-l 

SYNC. 

01 

LAST 

-JUL I AN-DAY 

USAGE 

COMP- 1 

SYNC. 

01 

STARTIM&-ST.O-OERTH 

USAGE 

C0MP- 1 

SYNC. 

01 

KEY- 

ERROR 

PIC 

S9999 

USAGE 

CO^P 

SYNC. 

«• 

These  are  used  when 

READER 

IS  called. 

01 

LINK 

-LATITUDF 

USAGE 

CO^P-l 

SYNC. 

01 

LINK 

-LONGITUDE 

USAGE 

COMP-l 

SYNC. 

01 

LINK 

-DYNAMIC-HEIGHT 

USAGE 

COMP-l 

SYNC. 

01 

LINK 

-eof-elag 

PIC 

S9999 

USAGE 

CO“P 

SYNC. 

*5ETK  E  Y  ROUTINE 

PROCEDURE  OIVISION  USING  F I RS T- JUL I AN-D A Y ,  L AST- JUL I AN-DA Y . 

starting-std-depth.  key-error. 

MOVE  ZERO  to  EOF-^LAG.  key-ERROP. 

OPEN  INPUT  OYNAMIC-HE IGHT-FILE. 

MOVE  STARTING-STD-DEPTH  to  NOM-STD-DEPTh. 
move  F IRST-JUL I AN-OA Y  TO  NOM-F I RST- JUL I AN-DA Y . 

START  OYNAMIC-HFIGHT-F ILE  USING  KEY  PECOPD-KEY  =  GENE  R I C-KF Y 
INVALID  KFY  MOVE  1  TO  KEY-EPROR 

CLOSE  OYNAMIC-HEIGHT-FILE. 

GO  TO  SETKEY-EXIT. 

SETKEY-EXIT. 

EXIT  PROGRAM. 


READER-ROUTINE. 

ENTRY  'READER*  USING  LINK-LATITUDE.  L I NK -LONG I TUDF . 

LINK-DYNAMIC-HE IGHT.  LINK-EOF-FLAG. 
MOVE  EOF-ELAG  TO  L I NK-EOE-FL A G . 


READ  DYNAMIC-HEIGHT-FILE 

AT  END  PERC OHM  CLOSE-UP 

GO  TO  EXIT-READER. 

if  standard-depth  >  STARTING-STD-DEPTH 
THEN  PERFORM  CLOSE-UP 
GO  TO  EXIT-READER 

ELSE  IE  JULIAN-DAY  >  LAST-JUL I AN-DAY 
THEN  PERFORM  CLOSE-Ud 

SURTRACT  90  from  latitude,  ISO  FROM  LONGITUDE  to 
change  THE  RANGE  from  0-180  AND  O-DGO  TO  -90  TO  *90 
AN->  -180  TO  *18  0. 

ELSE  SURTRACT  90  FROM  LATITUDE  GIVING  L INK-LAT I TUDF 

SURTRACT  180  FROM  LONGITUDE  GIVING  LINK-LONGITUDE 
MOVE  DYNAMIC-HEIGHT  TO  L I NK -D YN AM  I C-HF I GH T . 

EXIT-READER. 

EXIT  PROGRAM. 

CLOSE-UP. 

MOVE  1  TO  EOF-FLAG,  L I NK-EOF-FL AG . 

CLOSE  DYNAMIC-HEIGHT-FILE  . 


C* 

C* 

c* 

C* 

c* 

C4 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


.  PROGRAM  OYF'IT  * 

o 

WRITTEN  RY  KEN  MQONFY  * 

« 

CHANGES  * 

1)  Changed  FROM  main  routine  to  subroutine  to  rf  called  hy  A  * 

main  COBOL  ROUTINE  IN  ORDER  TO  FACILITATE  FORTRAN  AND  COBOL  * 
SUBROUTINE  COMMUNICATION.  * 

2 )  ALL  CALLS  TO  SPL4  HAVE  BEEN  DELETED  AS  SPLA  has  hFFN  a 

DELETED  FROM  THE  D YE  I T  SYSTEM.  « 

ENTERED  H/10/19S0  BOB  STARK  * 

3)  ADDED  STANDARD  DEPTH  OF  7ERO  METERS.  ENTERED  7/1/19P0  R.  STARK 


SUBROUTINE  nyFORT 
COMMON  A  (  1  HR ) 

COMMON  XUAT A ( 3. 1 000 ) . YDAT A (  1 000 ) , WDATA  (  1000  > 
COMMON  WORK (16) 

DIMENSION  7  (is) 

DATA  7/0.  ,  20  • 4  30. .SO. ,  75. • 100. ♦ 125. . 150. *200. ,250, 


. 300. .400. ,500. , 


A  s  0  0  • ,700, 

.500. ,900. , 1000./ 

C* 

DPDP 

STANDARD  OEp  TH 

c* 

ALAT 

CENTER  LATITUDE 

c* 

ALON 

CFNTE»  longitude 

c* 

CTTmE 

CENTER  IN  TIME  (JULIAN  DAY) 

c* 

YRNG 

range  for  latitude 

c» 

XRNG 

range  for  longitude 

c* 

TRNG 

range  FOR  TIME 

c* 

IPRT 

IF  1  DATA  PULLED  off  FILE  BY  SPL3  WILL  HE  PRINTED 

10  READ ( 1 «900«END  =  40 ) DROP , AL A T . ALON . CT I  ME , YRN^ , XRNG , TRNG . IPRT 
900  FORMAT  <F  5. 0.2X.  1X.F6.?»1X.F7.2,1X,F4.0,1X,FS.?.1X,F6.2,1X, 

1F4.0, IX , I  1 ) 

IF (DPDP.EQ.RR99. )  GO  TO  40 

PRINT  903,DPDP.ALAT, YRNG, ALON. xRNG.CTIME,TRNG 
903  FORMAT(*l  STANDARD  DEPTH  =  * ,F5.0,7X« 'LATITUDE  =  »,F6.?,«  ♦  /-  • , 


r 

c 


A  F5.2.7X, ’LONGITUDE  = 

R  Fto.O.  •  ♦/-  »  .  F  4 . 0 ) 

:  this  section  checks 

DO  20  1=1 » IS 
IF (DPDP.EQ. 7 (I)>  GO  TO 
20  CONTINUE 
PRINT  R  0  2 

r02  FORMAT ( '  DEPTH 
DPDP  WASN'T 
UNTIL  WE  FIND 
25  READ (1,901, END  =  4  0 ) 

rO 1  E  ORM A  T ( F  S . 2 ) 

IE  (DPDP.EQ. 99. )  GO  TO 
GO  TO  25 


* . E  7 . 2 , *  ♦/-  ' .F6.2.6X, • JUL I  AN  DAY  =  ', 

TO  MAKE  SURE  THAT  nPDP  IS  A  STANDAPD  DEPTH 


30 


ENTERED  IS 
A  STANDARD 
A  99  CARD 
DPDP 


NOT  A  STANDARD  DEPTH. • ) 

DEPTH.  NOW  WE  READ  AND  DUMP  RECORDS 
WHICH  MARKS  THE  END  OF  THIS  SET  OF  POINTS 


10 


30  CALL  SPL3(OPDP,ALAT, ALON, CTIME, XRNG, YRNG. TPNG.IPRT) 
GO  TO  10 


n  r> 


1 


4  U  PE  TURN 
END 

SUhROOT  INF  SPl  3  (  OPOP*  ALAT  ♦  ALON.CT  IMF  *  XHNG  ,  y  PNG  ,  TP  Nr,  ,  IPwT  ) 

COMMON  A ( 1  hR ) 

COMMON  XOATA(3»1000) »YUATA(1000)  ,  W  D  A  T  A  (10  0  0) 

COMMON  WOPK(lS) 

DIMENSION  X  min  {  ?  )  «  XMAX  (  ?)  .NODES  (2)  «  NOE  R  I  V  <  ?>  ,  COFE  (14)  .  x  (  ?  ) 

OAT  A  NO  I  M/?/, NODES  <  1  )  /4/. NODES  <  2  )  /4/  *  X  TRAP/  1  .  /  «  MCE  /  1  4/  ,  Nwrk  /  \  1  / 
kEAL  MINLAT.  MINlON,  MAXLAT.  maxlON 
I NTtGt P*P  EOF,  KEYfwP 

C  This  SECTION  REWRITTEN  O-13-1P80  BY  COOP  HOP  STARK  TO  PPOVIOF 
C  COMPATAH IL  I  TY  wITh  NEW  COBOL  INDEXED  SEGUE  NT  I AL  RECORD  RFTrIEVAL 
C  SUhPMOGPAMs  ADDED  TO  ThE  OYFIT  SYSTEM. 

XMlN ( P) = AL AT -YRNG 
MINlAT  =  XMI N ( ?) 

XMAX ( 7 ) =AL AT ♦ YPNG 
MAXLAT  =  X MAX(P) 

XM IN ( 1  )  =ALON- X  PNG 
MINlON  =  XM  TN  l  1  ) 

XMAX ( 1 ) =ALON»XPNO 
MAXLON  =  XMAa  (  1  ) 

TMIN=rT I M E-TP NO 
TMAX=CTIME*TPNG 
TmINP  =  0. 

Tmax?  =  0. 

NDAT A=0 

C  CHECK  FOP  JANUARY  TO  DECEMBER  WRAPAROUND  IN  START  oaTE. 

IF  (  TMIN.oE . 1)  30  TO  m 
TmIn=Tmin*36h. 

Tm I N?= 1  . 

Tmax?  =  Tma  X 
TMAX=  3GS. 

GO  TO  HI 

C  CHFCK  FOR  DECEMBER  TO  JANUARY  WRAPAPOUND  IN  The  FND  DATF. 
h  IF  ( Tmax.LF  .3GS.  )  GO  TO  in 
T  m I n  P  =  1. 

Tmax?  =  TMAX-3GS. 

TMAX  =  <  4  S  • 

C 

C  SET  ThF  STAPTTNG  FILE  xEY  TO  DESIRED  DYNAMIC  DEPTH 
C  AND  JULIAN  Day  (TMIN). 

10  CALL  SET*EY(Tmin,  TMAXi  DROP*  kEYERR) 

IF  (KEyEPP.NF . 1 )  00  TO  PO 

C  IF  WE  GOT  HERE »  ThE  SETKEY  ROIITINF  COULD  NOT  FIND  A  KEY  AT 

C  THE  STARTING  <EY  POINT.  NOW  wF  TRY  INCREMENTING  THE  START  DAY  WITH 
C  THE  HOPE  OF  FINDING  A  OOOD  KEY.  IF  wl  CAN'T  FIND  A  GOOD  KEY  FROM 

C  TMIN  TO  TMAX,  THF  AREA  MUST  HE  OUTSIDE  THE  RANGE  OE  T Ht  FILE. 

TMIN  =  TM I n* 1 . 

IF  (TMIN. LT. Tmax )  GO  TO  10 

0.  K •  wE  BOMHEO  OUT  ON  THAT  ATTEMPT  TO  SET  THE  KEY. 

CHECK  FOR  A  P  ENTRY  POINT  CONDITION  (IE  DECt *HER- JANUARY  WRAPAROUND) 


n  n  n  n  on  on  noon  r>  n  r>  o  n 


C  IF  F  OiJ'mD  ,  MHAMCii  TO  THE  SECONO  *F  Y  SETUP  ROUTINE 
IF  (T«IN2.NF.O.  .(jh.TMAX2.NI  .b.  )  GO  TO  2? 
r 

C  SOkky  -  CAM*T  SET  A  EILF  k  f  y  t 
PPINT  1S.DP0P 

IS  FORMAT  (t  Nr,  OATA  IN  FIlF  FOP  STD  DEPTH  =  ».Eh.O) 

C  WEAL)  AND  Dump  POINTS  UNTIL  a  00  IS  FOUND. 

1  7  PL  AO  030 . < (?) 

IF ( X ( 2 ) .NF .09. )  GO  TO  17 
GO  TO  70 
C 
r 

C  THIS  IS  the  MAIN  PEAO  LOOP. 

20  CALL  P  E  A  U  E  p  (ALATT,  along,  DYNHGT.  EOF) 

IF  (EOF.NE.1)  GO  TO  2S 

file  is  at  end 

CHECK  FOP  2  STAHT  POINTS. 

IF  (TWIN?. FO.n.. AND. TMAXP.EQ.O.)  GO  TO  40 

FlLt  DOES  HAVE  ?  START  POINTS  DUE  TO  DFCEM°EP- JANUARY 
WPAPAPOIJNO.  NOW  wF  have  TO  SET  ThF  KEY  FOP  ThF  SECOwO 
STAhT  POINT  OF  Thf  file. 

23  Thin  =  T  M  I N  n 
T  h  A  X  =  Thaxs 
TMIN2=o, 

T  M  A  X  2  =  0  . 
i,0  TO  10 

this  section  checks  the  data  fop  validity. 

THE  NEXT  STATEMENT  ADJUSTS  THE  LONGITUDE  IF  ThE  CENTEP 
IS  i-  aST  OF  THE  INTEPNAT  IONAL  DATE  LINt  AND  THE  POINT  IS  WEST. 
2S  IE  (MINLON.LT.-1HO. . ANO. ALONG. GT.O.)  ALONG  =-360. ♦ALONG 

IS  THE  NEXT  STATEMENT  ADJUSTS  THE  LONGITUDE  IF  THE  CENTEP  IS 
WEST  OF  THE  INTERNATIONAL  DATE  LINE  AND  ThE  POINT  IS  FAST. 

IE  (MA XLON.G r . 1  HO . .ANO. ALONG. LT. 0. )  ALONG  =  3h0.*AL0NG 

NOW  CHECK  THAT  LAT  AND  LONG  APE  pIThIN  RANGE. 

IF  ( ALA  TT.GT.MA  XL AT. OP. ALONG. GT.MA XlON.UP.ALATT.lt. MINE  AT. OP. 
1  ALONG.LT .mi \l ON)  GO  TO  ?0 

0.  K.  wE  MADE  IT.  DATA  IS  GOOD. 

G  DATA  =  N I)  A  T  A  ♦  ] 
vvOATA(NDATA)  =  1  . 

XlJATA  (2.NDAT4  )  =  ALATT 

X  D  A  T  A ( 1 , NO A  T  A )  =  ALONG 

YOATA(NUATA)  =  OyNHGT 
IF  (  Nr'  A  T  A  .  GE  .  1  9  0  0  )  GO  TO  4  0 
GO  TO  2  0 

RESUME  Oh.  KEn  moONEY’S  OWIGINAL  SPLS  POOTINF. 

40  IF(IPpt.NE.I)  GO  TO  So 
PP I  NT  P 1 0 

PlO  FOPMAT  ( 1H-. 1  OX . 3HLAT ,8X .4HL0NG.  SX.hHDY  HOT,//) 

PRINT  u?0  ,  (  I  ,  XDATA  (2,  I  )  ,  XOATA  (1  .  I  )  ,  YDAT  A  (  I  )  ,  1  =  1  ,NOATA  ) 

P20  FOPMATIIh  . I 3,SX ,E6.2,SX,F 7.2.  SX.Eh.3) 


r>  o 


50  C4lL  SPLCw  (NOIM.nDATA.XMTN.XMAX,  NODE  S.XTRAP.COEF  »  NCF  ,  N«i  Rk,  ,  I  ERROR) 
IF ( IERROR.GF.3?)  GO  TO  bO 
PRINT  940 

94  0  FORMAT ( 1 H 1 , ?X . 3HL A T *  8 X « 4HL0NG «  7X « .3RD A  Y  * b X  , 3HDIR.  7X  * 3hSPD, 7X  » 

1 1HU.9X, IHV#bX«5HN0ATA) 

READ  THE  LATITUDE  AND  LONGITUDE  AT  WHICH  THF  SURFaCF  CURRENT 
IS  TO  hE  COMPUTED. 

60  READ  9 30. X  (?)  ,X  (  I) 

930  FORMAT < Fb . ? . 1 X . F 7 .2 > 

IF (X (2) .EQ.99. )  GO  TO  70 
IF ( IERROh.GF.3?)  GO  TO  60 
C0N1=-1 .2/SINIX (2) /57.29b) 

C0N2  =  ?.4/SIM(2.t>X (?) /S7.296) 

NDERI V ( 1 ) =0 
NDERI V  (?)  =1 

VX=CON1«SPlOE (NDIM.X.NOERIV.COEF , XMIN.XMAX.NODES. I  ERROR) 

IF ( IERrOr.GF . 3?)  GO  TO  70 
NDERI V ( 1 ) =1 
NOERI V (2) =0 

VY  =  CON?*SPLD£ (NDIM, X ,NDER  I  V.COEF ,  XM  I  N  .  XM  A  X  .  NODES  .  I  ERROR) 

IF ( IERROR.GF. 3?)  GO  TO  70 
VEL=(VX*VX*VY*VY)*<MI.S 
ANG=57.296»ATAN? ( VX . VY) 

IF ( ANG.LT.O. ) ANG  =  ANG+  360 . 

PRINT  950.X (?)  ,X  (  I)  .CTIME, ANG.VEL, VX, VY.NDATA 
950  FORMAT ( 1H0.P6.2.5X.F 7.2 .5X .F 4.0 »5X.F 4.0 . 3(5X,F5.?).5x.T3) 

GO  TO  60 
70  RETURN 
END 

SUBROUTINE  SPLCW (NDIM.NDATA.XMIN.XMAX, 

1  NODES, XTRAP, COE F.NCF.NWR^ . IFRROR) 

C  DESCRIPTION  OF  SPLCW  IS  AVAILABLE  IN  NCAR  SOFTWARE  SUPPORT 
C  LIBRARY,  NCAR  TECHNICAL  NOTE  IA-105,  MARCH  1975. 


1. 1ST  OF  CAPTIONS 


Figure  I.  Node  arrangement  us<‘<l  in  .spline  fit  analysis. 

Figure  2.  Average  geostrophie  enrrents  in  the  Southern  California  Bight  for  the  month  of 
Jarman . 

Figure  3.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
Fehruarv. 

Figure  4.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
March. 

Fillin'  5.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
April. 

Figure  6.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
May. 

Figure  7.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
Ji  rie. 

Figure  8.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
July. 

Figure  9.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
August. 

Figure  10.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
September. 

Figure  11.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
October. 

Figure  12.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
November. 

Figure  13.  Average  geostrophie  currents  in  the  Southern  California  Bight  for  the  month  of 
December. 
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